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1 Introduction 

In [3], this author introduced the regularized sampling method. A method which is based on Shannon's 
sampling theory but applied to regularized functions, hence avoiding any (multiple) integration(s) 
and keeping the number of terms in the cardinal series manageable. It has been demonstrated that 
the method is capable of delivering higher order estimates of the eigenvalues at a very low cost. The 
purpose in this paper is to extend the domain of application of this method to the problem at hand. 

2 Main results 

Consider the following Sturm-Liouville problem with discontinuity inside a finite interval, 

-y" + q(x)y = /j, 2 y , x £ (0, vr) 

j,'( ) = = y(7r) , , 

y{d + 0) = ay{d - 0) 1 1 

y '(d + 0) = a~ V(d - 0) 

where a > 0, a ^ 1, < d < n and q G L 2 (0, n). 

It has been shown in [lj that the problem has a point spectrum and each eigenvalue has multiplicity 
one and accumulate only at +oo. The purpose in this paper is to compute the eigenvalues of ( 12.11) 
with the minimum of effort and a greater precision using the newly introduced regularized sampling 
method [3], an improvement on the method based on sampling theory introduced in [2]. 

Let y^x,^) and yn(x,fi), be the solutions of the base problems 

-y'L + Q( x )yi = n 2 yL , x e [o, d] , . 

y L (0) = 1, y' L (0) = [Z - Z) 
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and 



respectively. 



-y'k + q(x)vr = 1± 2 Vr 



x G [d, 7r] 



Vr{^) = , y' R (n) = 1 



(2.3) 



(2.4) 



Using y'iji, ji) = a, an unknown constant to be determined along the eigenvalue parameter fj 2 
the discontinuity conditions give 

ay R (d,/ji) = ay L (d,/ji) 
ay' R (d,fi) = a- l y' L (d,n) 

Note that a ^ since otherwise, yi(d,fi) = y' L (d, fi) = leading to |/l(x, /z) = which contradicts 
yi{0,fi) = 1. A necessary and sufficient condition for non trivial solutions is A(fi) = where the 
characteristic function A is defined by, 



AM 



ay Lid, n) yB.(d,n) 

a~ l y L {d,n) y R (d,v) 



(2.5) 



Thus the eigenvalues of the problem at hand are seen as the square of the zeroes of A. Let /j, k be 
such a zero, a will take the value 

ay L (d, n k ) 



yR(d,/J>k ) 

The eigenfunction associated to the simple eigenvalue Hk is defined by 



y[k]{x,Hk ) = 



yL{x,Hk ) if < x < d 
QikyR(x,fi k ) if d < x < n 



(2.6) 



(2.7) 



In the following we shall present some important properties of the functions yi(x, ji) and y R (x,/i) . 
We shall need the following well known lemmata, 

Lemma 2.1 sin z/z and cosz are entire as functions of z and satisfy the estimates 

\wnz/z\ < 7 e |7m zl /(l + \z\) and \cosz\ < e |/m ' A 

where 7 = 1.72. 



Lemma 2.2 If (3 is a positive constant and if is a positive function satisfying, 

pit 

z(x) <P+ <p(t)z(t)dt 

J x 

then 

z(x)<(3exp< / ip(t)dt 



Theorem 2.3 yL(x,/i), y' L (x,fi), y£,(x, //) — cos fix, y' L (x,fi) +fisinfix G PW X as functions of fi for 
each fixed x G (0, d] and satisfy the growth conditions, 



f x 

y' L (x, n) + /ism/ix — / cos/x(x — t)q(t) cos(/j,t)dt 

Jo 



\y' L (x, fi) +/isin/ix| < 7 4 e x|ImMl 



7l e*IWI 

72 


1 + 

7 


/i 7T 

i 


1 + 

7 4 e x 


h\tt 

Imfi 
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y R (x,fi), y' R (x, /i) + sm ^ x \ y' R (x,ii) — cos ^(tt — x) G PW W - X as functions of fi for each fixed 
x G [d, 7r) and satisfy the growth conditions, 

\y R (x,n)\ < — 2& e (*-*)l*H 



y R {x,n) + 



sin /i(rr — x) 



// 



< 76 (7r-x)|Im^| 



/,) - cos ^(tt - x) I < —iL—et*-* ^ 
where 71, • • -,77 are some positive constants. 



Proof: For t/iwe have, 



y L (x,fi) = cos fix + 



f 

Jo 



sin /i(x — t) 



Q(t)y L (t, n)dt 



and 



y' L (x : n) — —/ism /jlx + / cos fi(x — t)q{t)yi,(t, fi)dt. 

Jo 



(2.8) 



(2.9) 



Standard arguments show that yi, y'i are entire functions of \x for each x G (0, d\ and belong to PW X 
for each x G (0, d] as function of /i. 

\y L (x,Li)\ < e x ^ + r^ne^-W^ \q(t)\ ■ \y L (t,fi)\dt 

Jo 

e - x \im,\ \y L{Xjfi)l < 1+ [ X lo n \q(t)\ ■ e-^ \y L (t,fi)\dt 

Jo 

using Gronwall's lemma gives, 

e-«|^(x,/i)| < exp| 7 ovr^\(t)|^ 



< exp ^ 7 7T ^ k/( /")| (It 



PIT 

^0 



so that, 



\Vl(x,h)\ < exp {to* jT \q(t)\dt} e x ^ = 1± e x ^ 
where 71 = exp {^ q tx \q(t)\dt}. 



(2.10) 



s sin fji(x — t J 
y L (x,fi) - cos fix = I q(t)y L (t, fj)dt 

t 1 



(2.11) 
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\ul(x, jj) — cos fxx\ < 



sin /i(x — t) 



< 



< 



/i 

sin n(x — t) 



7o 



H(x-t) 
(x-t) 



1 + \/jl\(x — t) 



\q(t)\-\y L (t,fi)\dt 
x-t)- \q(t)\-^ Im ^dt 
q{t)\- lie tlIm ^dt 



7T 



1 + |/i|7T 



(2.12) 
(2.13) 
(2.14) 



where 72 = 7T7o7i J* \q(t)\dt. 
Now, 



/•x px 

y' L (x : /j,) + /ism/ix — / cos /jl(x — t)q(t) cos(iit)dt — / cos fi(x — t)q(t) {yi(t, fi) — cos fit} dt (2.15) 
Jo ./o 

so that, 

rx 

y' L (x, fx) + \x sin/ix — / cos /j,(x — t)q(t) cos(/j,t)dt 

Jo 

< [ e^-W^ \q{t)\ \y L (t,n) - cos fit\dt 
Jo 

\q(t)\dt 



1 + |/x| 7T 

where 73 = 72 /* |g(f)| dt. 

/i) + /i sin /zx| 



< 



73 e x|Im/a| 

1 + \n\ir 



cos /z(x — t)q(t)yL,(t, /i)dt 



P7T 

lie x\i*w\ / | g (t)|dt = 74e *IWI 
Jo 



(2.16) 



(2.17) 
(2.18) 



where 74 = 71 |g(t)| dt. 
As for we have, 



sin/i(7r — x) f w sin fi(x — t) 



Vr(x,h) = -^ITXl rz _ / — r-v- - q (t)y R (t,n)dt 

fi J x fl 



(2.19) 



and standard arguments show that y# and are entire functions of /i for each x G [d, n) and belong 
to PWn-x for each x G [d, ir) as function of /i. Also, 



Vr{x,h) + 



sin /x(7r — x) Z" 7 " sin /i(x — i) 



// 



J X 



q(t)y R (t, n)dt 



(2.20) 



and 



y' R (x,fi) - cos jj,(Tt - x) = / cos n(x - t)q(t)y R (t, n)dt (2.21) 

J X 

are entire functions of \i for each x G [d, 7r) and belong to PW^-x for each x G [d, 7r) as function of \x. 



Using the Lemma 2.1. (12.191) gives, 



e -(^)|im M | < 7o ^ + / 70 7r|g(t)| ■ e-^^\y R (t, n)\dt (2.23) 

Gronwall's inequality (Lemma 2.2) yields, 



e-^ Im ^\y R (x, f i)\< lo -——exp{ lo n / \q(t)\dt\ (2.24) 

from which we get, 

\y R (x,fi)\ < 7o ^_exp( 7o vr r|g(t)|^le^l Im ^ < e (— )l^l ( 2 .25) 

where 75 = 7 vrexp {7 tt \q(t)\dt} . 

sin/i(7r — x) f* sin u(x — t) ,. , . , 

x, /i) + — ^ L = - / — ^ L q(t)y R (t, fj,)dt 

H J x fi 



sin/i(7r — x) 
y R {x,n) + - 



where 76 = 70^75 J™ \q{t)\dt. Similarly, 



" \ [0 l + |/i|7rl + |/i|7r7 m;i J + 

(2.26) 



|^(x, M )-cos/x(7r-x)| < re^l Im ^|g(t)|-|^(t,/i)|rft 

£ {Trfcr i?(i)l *} e< ^ ,llwl 

= ^ (w-a.)|ImM| C2.27) 

l + |/i|vr V ; 

where 77 = 75/* |g(t)|dfc. 

Thus, y R (x,fj,), y R (x,fj,) + sm v(n-x) ^ y' R ^ x ^^ _ cos ^ i (7 r _ a;) are entire functions of /x for each 
x G [d, 7r) and belong to PW^-x for each x G [d, tt) as function of /x. 

Although we have obtained much higher estimates in [5] and jl] to the expense of subtracting 
terms involving multiple integrals, and as in [3J, we shall stick with the estimates given in Theorem 
2.2, hence avoiding any (multiple) integration(s) and show by the same token that we can get a 
higher order estimate of the eigenvalues of the problem at hand at a very low cost. In fact we do not 
have even to keep on increasing the number of sampling points. 

Let PW a denote the Paley- Wiener space 



PW a = {/entire, |/(/i)| < Ce^ lm ^, [ 



|/(/i)| 2 d/i < 00} 

R 
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Let hki be defined by 



h n (fx) = 

h 2 (fi) = 

h 21 {fi) = 

h 22 (/-0 = 



sin 6fi 
8/1 

sin 9fi 
9fj, 



(y L (d,fi) - COS fid) 
(y' L (d,fi) +fi sin fid) 

nf 

0/1 



^ (y' R (d,fi)- cos fi(n-d)) 



Then we rewrite y L (d,fi), y' L (d,fi), yii(d,fi) and y' R (d, fi) as 

y L (d,fi) = h n (pi)[^) +cosfid 
y' L (d,fi) = h 12 (fi)(^y m -fismfid 
y R (d,fi) = h 21 (fi)[^y m -EE^pS 

y' R (d,fl) = h 22 (fl) {^y" 1 + COS fl(7T-d) 

Theorem 2.4 Let d be a positive constant and m be a positive integer (m > 2). The functions 
hki, (k,l = 1,2) belong to the Paley space PW a with a = cr + m6 and satisfy the estimates 



IMaOI < 



7 



n\lmu 



(l + %l) m 

k, I — 1, 2 /or some positive constant 7, <x = max {d, vr — d}. 

Proof: That h^i are entire and satisfy the given estimates is a direct consequence of Theorem 
2.2 and the fact that ^jj^ is an entire function of fi and satisfy the estimate in Lemma 2.1. 

Since the h k i(fi) belong to the Paley- Wiener space PW a for each k, I — 1, 2, they can be recovered 
from their values at the points fij = j^, j <E Z , using the following celebrated theorem, 

Theorem 2.5 (Whitaker-Shannon-Kotel'nikov) Let h E PW a , then 

^ sin a(fi - fij) 

//j = j^. T/ie series converges absolutely and uniformly on compact subsets of C and in L 2 d ^R). 

For all practical purposes, we consider finite summations, therefore we need to approximate hki by 
a truncated series The following lemma gives an estimate for the truncation error. 

11c 

h k i(fi). Then, for \fi\ < Nir/a, 



Lemma 2.6 (Truncation error) Let h^P(fi) = J2f=~N (l^j) s '°^1 - denote the truncation of 



< 



sva.fi\c 2 



7r(7r/o-) m -Vl - 4" m+1 
where c 3 = ||/i m_1 /iM(/i)||2- 



VW^J 



+ 



a)- fi 



a) + fi 



(N + l) 



m—l 
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Proof: Since /i m_1 /i£;/(/-0 £ L 2 (— 00,00), Jagerman's result (see [TT], Theorem 3.21, p. 90) is 
applicable and yields the given estimate for the hki, k, I — 1, 2. 

An approximation Bn to the characteristic function B is provided by replacing the hki by its 
approximation h^p, and we obtain at once, 

Lemma 2.7 The approximate characteristic function B N satisfies the estimate, 



\B(v)-B N (n)\ < 



sin 9fi 



9a 



sin /i|c 4 



7r(7r/o-) m -Vl - 4" m+1 



+ 



y/(Nn/a) - n ^(Nrr/a) + li 



(JV+l)m-l 



/or some positive constant C4. 
We claim the following, 

Theorem 2.8 Let~p 2 be an exact eigenvalue of B (multiplicity 1) and denote by [i 2 N the correspond- 
ing approximation of a square of a zero of B N . Then, for \[In\ < Nn/a, we have, 



\/j.n - v\ < 



\B'tf)\ 



sin 9fi n 



9fi 



N 



sin /i7v|c 4 



+ 



7r(7r/cr) m -Vl - 4- m+1 
1 



^(Nn/a) - fi N y/(Nn/(r) + li n 



(N + l)( m_1 ) 



where the inf is taken over a ball centered at lin with radius |///v — lA a> n d not containing a multiple 
ofn/9. 

Proof Since ~fi is a simple zero of B , we have, 

B(p)-B(fi N ) = (-p-fM N )B'(Ji) 

for some fx. Thus, 

\B(p) - B{» N )\ 



If 1 ~ Mjv I 



< 



|s'0Z)| 

sin 9 ll jy 



|£'(pOI 



1 



Sin ;Ujv|C4 



7T 



+ 



(7r/cr) m -Vl - 4- m+1 



y/(N7T/a) - fi N ^{Ntx/o-) + n N 



(N+r 



m—1 



where the inf is taken over a ball centered at /i^ with radius \hn — lA and not containing a multiple 
of n/9. Thus, the result. 
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3 Numerical examples 



In this section, we shall work out an example to illustrate our method. We shall take N = 40 , 
m — 6. We have taken 9 = cr /(N — m) in order to avoid the first singularity of ( s 'g^ JV j ■ The 

sampling values were obtained using the Fehlberg 4-5 order Runge-Kutta method. 
Example 3.1 

-y" + q(x)y = /j, 2 y , x G (0, vr) 

y'(0)=0 = y( 7T ) 

y{d + 0) = ay{d - 0) y ' 

y'{d + 0) = or V(d - 0) 

where a = 2, d = 1 and = x. 



Index 


Exact 


RSM 


Absolute Error 


Relative Error 


1 
2 
3 
4 


1.22788546912 
1.83749384727 
2.68396812434 
3.85661744715 


1.227885469249 
1.837493847255 
2.683968124476 
3.856617447367 


1.31357098 x lO" 10 
1.6491398678 x lO" 11 
1.32813906505 x lO" 10 
2.1294955259 x lO" 10 


1.06978298879 x lO" 10 
8.9749408974 x 10~ 12 
4.94841594058 x 10~ n 
5.52166647358 x lO" 11 



The graphs of the first four exact eigenf unctions and their approximations are displayed in Figure 
1. The inner products of any two different approximate eigenf unctions are of the order of 10~ 10 . 



4 Conclusion 

In this paper, we have used the regularized sampling method introduced recently [3] to compute the 
eigenvalues of Sturm-Liouville problems with discontinuity conditions inside a finite interval. We 
recall that this method constitutes an improvement upon the method based on Shannon sampling 
theory introduced in [2] since it uses a regularization avoiding any multiple integration. The method 
allows us to get higher order estimates of the eigenvalues at a very low cost. We have presented an 
example to illustrate the method and compared the computed eigenvalues with the exact ones. We 
have also computed the approximate eigenfunctions and compared them with the exact eigenfunctions 
and checked that the inner products of any two different approximate eigenfunctions are of the order 
of 10 -10 , a point which validates our method. We shall present in a future paper a generalization of 
the above result together with extensive numerical computations. 
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